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Abstract 

Using an analytical model for the string network we show that the 
kurtosis of cosmic microwave background (CMB) temperature gradi- 
ent maps is a good statistic to distinguish between the cosmic string 
model and inflationary models of structure formation. The difference 
between the stringy and inflationary value for the kurtosis is inversely 
proportional to the angular resolution and to the number of strings per 
Hubble volume of the strings' scaling solution. If strings are indeed re- 
sponsible for CMB anisotropics then experiments with resolutions of a 
couple of arcminutes or smaller could determine it using this statistic. 
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1 Introduction 



Presently there are two main models for the formation of large-scale struc- 
ture. In the one, quantum fluctuations originating in an inflationary epoch 
become classical density perturbations (Guth & Pi 1982; Hawking 1982; 
Starobinsky 1982; Bardeen, Steinhardt & Turner 1983), which grow by grav- 
itational instability into the structures seen today. In the other, topologi- 
cal defects arising in a phase transition in the early Universe act as seeds 
for galaxy formation (Kibble 1976; Vilenkin 1985). Galaxy redshift sur- 
veys, measurements of peculiar velocities of galaxies and of temperature 
anisotropies in the cosmic microwave background (CMB) are being used to 
test these theories. But so far experimental data have not been sufficient to 
rule in favour of one model or the other, being in more or less good agreement 
with both (see e.g. Brandenberger 1992). 

Our goal is to find, using an analytical approach, a signature of cosmic 
strings on the CMB which is unique to them and not predicted in inflationary 
models. Cosmic strings are one-dimensional topological defects which form 
when a symmetry is spontaneously broken to a smaller subgroup such that 
the first homotopy group of the vacuum manifold is nontrivial. They are 
very thin lines of trapped energy with mass per unit length /x; the spacetime 
around them is flat, but with a deficit angle of a — 87rG/i (Vilenkin 1981). 

Inflationary models predict a gaussian distribution for the CMB temper- 
ature anisotropies measured at each angular scale (Efstathiou 1989). Tem- 
perature anisotropies from cosmic strings have been shown numerically by 
Gott et al. (1990) and analytically by Perivolaropoulos (1993b) to be also 
nearly normally distributed, in spite of the inherently non-gaussian nature of 
the effect of a single string: it causes linear steplike discontinuities in the mi- 
crowave sky (Kaiser & Stebbins 1984) (for a study of the statistics of density 
perturbations due to point-like seeds see Scherrer & Bertschinger (1991)). 
By the central limit theorem the combined effect of all strings results in 
a gaussian signal for a large number of strings . However, the probability 
distribution for temperature gradients at small enough angular scales does 
show a departure from normality (Gott et al. 1990). We use a model for 
the string network to calculate the gradient probability distribution (actually 
the moment generating function which is equivalent to it) and show that its 
kurtosis is a good statistic to discriminate between strings and inflation in 
CMB measurements at angular scales below a few arcminutes. The model 
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is also implemented as a Monte-Carlo simulation, which is used to obtain 
actual CMB maps from cosmic strings and check the analytical results. 

The reason why the distribution function for temperature gradients is 
more nongaussian than that for the anisotropics is that a superposition of S- 
function perturbations, the gradients of step-functions, approaches the gaus- 
sian distribution much more slowly than the superposition of the step-like 
temperature perturbations themselves. 

We are making use of an analytical model (see e.g. Perivolaropoulos 
(1993a,b)) for the temperature anisotropies caused by strings which the 
photons encounter between the time of last scattering and today. This 
model was applied to the study of peculiar velocities by Vachaspati (1992) 
and Perivolaropoulos & Vachaspati (1993)). It has been used to calculate 
the power spectrum and the amplitude of CMB anisotropies, as well as 
their probability distribution and corresponding moments. The temperature 
fluctuations were found (Perivolaropoulos 1993a) (see also Hara, Mahonen 
and Miyoshi (1993)) to be well approximated by a scale- invariant Harrison- 
Zel'dovich spectrum normalized to a value consistent with observations of 
large-scale structure (Perivolaropoulos, Brandenberger & Stebbins 1990) (for 
other studies of large scale structure formation by cosmic strings see Hara 
& Miyoshi (1993); Stebbins et. al. (1987); Vollick (1992)). The distribution 
function of temperature anisotropies was found (Perivolaropoulos 1993b) to 
be approximated by a gaussian to an accuracy of better than one percent. 
These studies were useful tests of the consistency of the cosmic string model 
with existing CMB data, but did not provide a way to distinguish it from 
other consistent theories based on inflation. 

An alternative to our approach is a full numerical simulation of the evo- 
lution of the string network, and a subsequent numerical evolution of the 
photons through it (Bennett, Stebbins & Bouchet 1992). An advantage of 
our analytical approach is that the explicit dependence on string parameters 
and angular scale can be shown, and the statistical fluctuations around the 
predicted values for actual experiments, due to having only a finite number 
of measurements, can be calculated. 
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2 Theoretical Model and Experimental Setup 



According to the scaling solution for cosmic strings (Bennett & Bouchet 
1988; Allen & Shellard 1990), there is a fixed number 10 of strings with 
curvature radius of the order of the horizon per Hubble volume at any given 
time. Their orientations and velocities are uncorrelated over distances larger 
than the horizon, and they move with an rms velocity of v s = 0.15c. In our 
model the strings are approximated as straight, so that the Kaiser-Stebbins 
formula for the temperature anisotropy due to the deficit angle of spacetime 
around the string can be used (Stebbins 1988): 

5T/T = ±47rGfx | k ■ itf a v a x e s ) | (1) 

where k is the direction of observation, v s is the velocity of the string, e s 
its orientation, and 7 S = (1 — v s /c)~ l l 2 . Photons passing in front of the 
string are redshifted, while those passing behind it are blueshifted. The 
deficit angle extends only out to a distance of one Hubble radius (H^ 1 ) from 
the string; this is due to compensation (Traschen, Turok & Brandenberger 
1986; Veeraraghavan & Stebbins 1990; Magueijo 1992): background matter 
redistributes itself when strings are formed so that the monopole and dipole 
moments of the energy density prturbations are zero. Therefore photons 
passing a string at a distance larger than the horizon scale are not affected 
by it. This cutoff is quite sharp because of gravitational shock waves caused 
by the compensating underdensities (Magueijo 1992). 

The time between now and last scattering is divided into N Hubble times 
ti with t i+ i = 2tj, and N = log 2 (to/ti s ) m 15 for a redshift of last scattering 
zi s = 1000. The apparent angular size of a Hubble volume at time ti is 

]/2 1/3 

given by Qm ~ \ ~ V f° r large redshifts in the matter dominated era 
and assuming Q = 1. Therefore, Hi+1 = 2 1 ^ 3 6 Hi , and 9 Hls = z^ 2 rad = 
1.8°. The strings encountered by the photon at subsequent Hubble times are 
assumed to be uncorrelated. This is justified because after a time t« a photon 
has traversed a large part of its original Hubble volume, so that at 2tj its 
new Hubble volume is mostly new space. Moreover a string moves a distance 
v s ti = O.lBcti which is approximately equal to the strings' mean separation, 
so that it is likely to encounter another string and intercommute with it. 
At each ti we assume the M string segments (each of length 2H^ X since the 
strings' orientations and velocities are uncorrelated over larger distances) per 
Hubble volume to have random positions, orientations and velocities. 
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The effect of all strings together is taken to be a superposition of the 
effect of the individual long strings. Finally, the effect of loops is considered 
unimportant compared to that of long strings, and the initial temperature 
inhomogeneities at the surface of last scattering are assumed to be negligible 
compared to those induced by the string network at later times. 

Neglecting the effect of loops can be justified as follows: The typical loop 
size is about 10~ 4 of the Hubble radius (Bennett & Bouchet 1988 ; Allen & 
Shellard 1990). At recombination this subtends an angle of 0.7". Only at a 
redshift z ~ 2 the angle becomes equal to 18", the smallest scale which we 
consider in this paper. Thus, the loops are too small for the Kaiser-Stebbins 
effect to apply. Instead, loops contribute to 5T/T via the usual Sachs- Wolfe 
integral (Sachs & Wolfe 1967), and the result of this integral is given by 
ST/T = where <3> denotes the relativistic potential at t\ s . Since on scales 
corresponding to a beam width of b> 18" many small loops contribute to 
$, the combined effect can (by the central limit theorem) be expressed as 
gaussian noise. We will estimate the magnitude of this noise in Section 6 and 
demonstrate that it is indeed negligible, in agreement with the conclusions 
of Stebbins (1993). 

Individual cosmic strings produce linelike discontinuities and plateaux in 
the microwave background. These are features which are most easily seen 
in temperature gradient maps. Experiments could either make line surveys 
or map out the temperature over a certain region of the sky, say a square of 
size 9° x 9°. Small angular resolutions b provide a better chance of seeing 
these features. We will quantify this statement in the following sections by 
calculating the moment generating function for the temperature difference 
between neighbouring areas that can be resolved by the experiment. 

Fig.l shows the result of a simulation for the temperature anisotropies 
expected in an experiment with b = 80" and 9 = 2.2° (i.e.lOOx 100 pixels). At 
each Hubble time tj rii string segments of length were placed at random 
over an area of size (9 + 9 H .) 2 (see Fig. 2). By the scaling solution 

m = M(0 + Hi ) 2 /P Hi (2) 

Note that the discontinuity lines seen in the map are not all at the position of 
the strings themselves. Each string at each of the 15 Hubble times gives rise 
to five such lines, one at its position, two at a distance H^ 1 parallel to it, and 
another two because it has finite length. Since strings are one-dimensional 
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objects, a square region with sides equal to one quarter of the horizon is 
crossed by about M/4 strings. By equation (1) the string shown changes the 
temperature of region A by (3r , and that of region B by —f3r, where 

(3 = AirGfjLjsVs (3) 

and r = \k • (v s x e s )\, v s — v s /v s - The direction of observation k is approxi- 
mately constant over the whole observed area, and e = e s x v s is a random 
unit vector since the strings' orientations and velocities are random, so that 
r = \k ■ e\ is uniformly distributed over the interval [0,1]. Each beam is 
assigned the temperature averaged over the area covered by it. 



3 Moment Generating Function for Temper- 
ature Gradients 

Let X be the random variable (RV) whose realizations are the possible values 
for the temperature differences between two neighbouring beams (pixels). 
Our goal is to derive its distribution function. We do so by calculating the 
moment generating function Mx(t) of X. The moment generating function of 
a RV is an important concept in statistics, it is defined by Mx{i) =< e >, 
where <> denotes the expectation value. The jth moment fij =< X j > of 
X is obtained from M x (t) by differentiation: 

ft= &L M * (t) (4) 

An inversion formula gives the (cumulative) distribution function Fx(x) = 
P(X < x), the probability that X < x (Lukacs 1960): 

F x (x + 5) - F x (x - 5) = lim - f *^e~ itx <\> x {t)dt (5) 

provided that x + 5 and x — 5 are continuity points of F. <f>x(t) = Mx(it) is 
called the characteristic function of X. Equation (5) is the general formula for 
a RV that need neither be absolutely continuous nor purely discrete. If it is 
continuous, it can be shown from (5) that a probability density p(x) = F'(x) 
exists and that it is the Fourier transform of <f>x(t)- 
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In the remaining part of this section we will derive Mx(t). The reader 
who is not interested in the details is advised to jump to equation (21) and 
continue reading from there. 

Since the temperature difference results from the superposition of the 
effects of the individual strings, we can decompose X into X = Y^\=\ X 1 , 
where X 1 is the RV for the temperature difference due to all the strings at 
Hubble time £«. In turn we can write X 1 as X 1 = YJjLi X}, where X 1 - is the 
RV for the temperature difference due to the jth string at time U. 

Now the moment generating function has the nice property that Mx+y (t) = 
Mx(t)My(t) if X and Y are independent RV. Since the Xj are independent 
in our model, Mx(t) can be written as 

15 ra» 

M x (t) = HM xi (t) , M xi (t) = Y[M xi (t) (6) 

i=l j=l 3 

tx i 

Therefore only M X i(t) =< e i > has to be calculated, which is done below 

by determining the possible values for XI- together with their probabilities. 

The two beams have the same temperature if they both lie outside the 
string's region of influence or both in region A or B (see Fig. 2). This is 
true because the observation direction k is almost exactly the same for both 
beams for b < 9 Hls . We are only interested in such small angular resolutions 
because the simulations have shown that for larger ones the sought-after 
non-gaussian signature disappears. Moreover we will find from the analytical 
result that the non-gaussian signature goes to zero as b tends to 9h 1s - The 
beams have unequal temperatures if either the string itself or one of the 
lines of discontinuity at Hubble volume cutoff (see dotted lines in Fig. 2) pass 
through either or both of the beams. 

First consider the string itself. The probability for it to cross the right 
beam (see Fig. 3 and Fig. 2) is 

p (V>H t + e)/i>) 2 ( ' 

for b 9 If it crosses the right beam it has uniform probability to be a 
distance n|,n G [0,1], away from its centre, i.e. to be tangent to a circle 
C n of radius n| around the centre. If it is tangent to C n to the right of the 
beam's centre and such that it does not cross the left beam (i.e. somewhere 
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along the arc between points A and B), then the possible values for Xj are 
(for r=l) 

X] = ±P(n-l) (8) 

each with probability p*y, where pi is the probability for a string which is 
tangent to C n to be so between A and B. To obtain Eq.(8) we have used that 
b <C Oh an d approximated the circular beams by squares: 



X) = ±l^(Kl - n-)(3 - b( b - + n-)(3) - {-(3)] = ±0(1 - n) (9) 



From Fig. 3 it can be seen that 

«i 1 . 1 — n. 

Pi = — = — arccos( — - — ) (10) 

7T 7T 2 

Similarly one obtains 

X; = ±f3(n+1) (11) 

if the string is tangent to C n to the left of the beam's centre such that it does 
not cross the left beam (i.e. between C and D), each with probability p 1 ^, 
where 

1 ( n + l \ <m\ 

p 2 = - arccos(— — ) (12) 

71 Z 

Finally, if the string is tangent to any other point of the circle, it passes 
through the left beam, and it has uniform probability to do so at distances 
m| away from that beam's centre, m G [0,1]. For a given m the possible 
values are 

X) = ±(3{n + m) and±/3(n-m) (13) 

each with probability p lE f, and p 3 = 1 — (p± + p 2 ). Since we could have 
started with the other beam, all the probabilities given above for the values 
of Xj have to be multiplied by 2. 

For each of the other four lines of discontinuity the same results apply, but 
with (3 replaced by (3/2. The probability for X 1 - = is 1 — 2^(1+4) = 1 — 10p\ 
Finally, r is uniformly distributed over the interval [0,1], so Xj — > rXJ. 

Putting all of this together (and setting f3 — 1) we arrive at 

M X i(t) = 

1 - lpp* + 2pm/ 1 dr f 1 dm f 1 d n ^(e (n+m)tr + 
jo Jo Jo 4 
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^—(n+m)tr _|_ ^(n—m)tr _|_ ^—{n—m)tr _|_ 

^ e (n+m)tr/2 _|_ e -(n+m)tr/2 _|_ e (n-m)tr/2 _|_ e -(n-m)tr/2^ _|_ 

+ e (-(n-l)tr + 



/V I' dn{ Pl( e (n-l)tr 

Jo Jo 2 

4( e («-l)W2 + e -(n-l)tr/2^ + 



P2( e (n+l)tr + e -(n+l)tr + 4 ^(n+l)tr/2 + e -(n+l)tr/2^|j 

We approximate £>i,£>2 and ^3 as constant over the interval < n < 1. For 
Pi = £>3 = 0.4 and P2 = 0.2 the integrals can be done with the result 

M X i(t) = 2p i {0A^{-smh 2 {t) - 16sinh 2 (^)) + 

i(0.6shi(2t) + 5shi(t) + 1.6shi(-))} + 1 - lOp* (15) 

where shi(t) = Jq dx sm ^ x . The next step is to calculate 

M xi (t) = fj M xi (t) = (M xi (t)) ni (16) 

since the Xj are identically distributed. Define 

A* = p'rii Ki M-^— (17) 

Hi 

Then 

Mxi (t) = (l + 2X V(t)-$) y _^ e 2A'(/(*)-5) (18) 

for n l large, which is the case here (f(t) is the function inside the curly 
brackets of Eq.(15)). This approximation makes it possible to do the final 
sum over Hubble times: 

M x (t) = J] M xi (t) = e 2 S!=i^(/(*)-5) (19) 
i=i 

Since 

15 15 1 Mb 



8 



the final expression for the moment generating function is 
M x (t) = 

e 9.4A(0.4^(-sinh 2 (t) - 16sinh 2 (|))) 
= 9.4A(±(0.6shi(2t) + 5shi(f) + 1.6shi(|)) - 5) ^ 
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4 Probability Distribution and Kurtosis 

Fig. 4 shows a plot of frequencies of temperature differences for b = 18" 
obtained from the moment generating function Mx(t) via the inversion for- 
mula Eq.(5) and from a simulation. The two curves are in good agreement. 
A gaussian distribution is shown for comparison. 

The kurtosis of the distribution of temperature gradients is a good statis- 
tic to distinguish between the string model and inflationary models. The 
kurtosis, defined as 

k ' ~ < (x- < x >y >2 {22} 

measures the peakedness of a distribution and how far out its tails extend. 
But this is how a gradient distribution caused by strings differs from a gaus- 
sian one (which is predicted by inflationary models): The sharp disconti- 
nuities along the location of the string produce a higher number of large 
gradients, while the plateaux give a higher number of small gradients than 
would be expected in the gaussian case. 
For a RV with zero mean 

h = ^ (23) 

where /ij is the jth moment given by Eq.(4). A Taylor expansion of Mx(t) 
around t — gives 

14 

k A = 3 + (24) 
A 

(see Fig. 5). The gaussian value for the kurtosis is 3. Recall that A = M-^—. 

"Is 

So the departure from the gaussian value is inversely proportional to angular 
resolution and to the parameter M describing the string network. This is the 
main result of our paper. For z\ s = 1000 and M = 10 we have 

4 ^ b (in arcminutes) ^ ^ 
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This shows that for beam sizes of a few arcminutes or smaller the string 
model predicts a non-gaussian kurtosis, whereas for larger beam sizes values 
very close to the gaussian one are predicted; consequently the latter exper- 
iments cannot be used to discriminate between strings and inflation using 
this statistic. Experiments with small angular resolutions have not yet de- 
tected CMB anisotropies. In the next section we discuss some of them to 
see whether they should in future measure a nongaussian value for A; 4 if the 
anisotropies are caused by strings. 



5 Confrontation with Observation 

Experiments only make a finite number of measurements and therefore do 
not measure the true moments. In order to distinguish between different 
models the distribution of the measured quantities around the exact ones 
has to be determined in each model. 

An experiment as described in Section 2 makes n = (|) 2 measurements 
of temperature differences between beams lying next to each other along a 
particular direction. We can interpret the ith measurement as one realization 
of the RV where is the RV X of the previous sections for the ith 
pair of beams. If we define 

y '- <g::w •^■-• n (26) 

then the are identically distributed, - Y^=i is the measured kurtosis 
kf, < yw > is the true kurtosis k± of the underlying distribution, and 

Var(Y^) =< (V«) 2 > - < > 2 = k 8 - k\ (27) 

where k$ = ^§ is the normalized eighth moment. By the central limit theorem 
the RV 

ri =1 (Y^-<Y^ >) 
y/ny/Var(YW) 

is normally distributed with unit variance for large n (Feller 1971), i.e. 

E « _ (yW_ < y(0 >) \ 1 re 

— — - < e = $(e) , $(e) = -= / dxe x 11 
^nyfVar(YW) J V2n 
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* P(kT~h <e) = *LEJ^\ (29) 



P 



From Eqs.(27) and (29) we can see that the measured kurtosis is normally 
distributed around the true value with variance a 2 = hs ki . For a gaussian 
distribution fcf = 105, while for a stringy distribution the Taylor expansion 
of M x (t) to eighth order gives 

29 1.8 0.017 
h = 105 + — + — + (30) 

Fig. 6 shows the distribution of the measured kurtosis about the predicted 
values in the two models for an experiment by Fomalont et al. (1990), who 
map out two areas each of size 7' x 7' with angular resolutions ranging from 
10" to 80". For a resolution of 18" a clear nongaussian signal is predicted, 
while for b = 80" the number of measurements is not sufficient to discriminate 
between the two models. If n were larger however, it would be possible to do 
so at 80". 

A second experiment we want to discuss is one proposed by the Cambridge 
group which is to start construction in the second half of 1993 (Lasenby 1992). 
It intends to map out an area of size 10° x 10° at resolutions ranging from 10' 
to 2°, and it was planned, among other things, to look for possible signatures 
from topological defects. Fig. 7 shows that for its highest resolution of 10' 
and M — 10 it would not be able to tell whether the underlying distribution 
is gaussian or due to strings. It would be nice if this new experiment could 
be designed also to make measurements at the scale of 1 arcminute, in which 
case it would be able to rule in favour of or against strings. Fig. 7 also shows 
a curve for a different value of the string scaling solution parameter M. If M 
should indeed be as low as 3 then strings would give a nongaussian kurtosis 
even at 10'. 

Based on our statistic, the Owens Valley experiment (Readhead et al. 
1989) with a beam width of 108" could detect a signal for cosmic strings if 
a map of contiguous patches were produced. Also, Melchiorri et al. (1993) 
are planning a 3 m telescope balloon experiment with a beam width of 2' 
which will have the potential of detecting non-gaussian signatures from cos- 
mic strings. 
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At 10' the nongaussian signature has all but disappeared, so that for 
experiments done at degree scale or at several degrees (e.g. COBE where b m 
7° ), strings do not predict a nongaussian distribution of either temperature 
gradients or temperature values. 



6 Effect of Noise 

Gaussian noise diminishes the nongaussian signature expected from strings 
alone. The moment generating function for the temperature difference be- 
tween two neighbouring beams due to gaussian noise is 

s 2 t 2 

M Y (t) = e~ (31) 

where s 2 has to be determined for each kind of noise. The stringy moment 
generating function with (3 = ATiGjj,^ s v s not set equal to 1 is 

M x (t) = 

e 9.4A(0.4^(-sinh 2 (/%) - 16sinh 2 (f ))) 
xe 9.4A(i(0.6shi(2/%) + 5shi(/ft) + 1.6shi(f )) - 5) ^ 

The measured temperature difference is given by X+Y, and since X and Y are 
independent, Mx+y{t) = Mx(t)My(t). Without noise the second moment 
of the probability distribution for temperature differences is /x 2 = 5.43A/3 2 
and the kurtosis is A; 4 = 3 + With gaussian noise the second moment is 

/4 oisc = 5.43A/3 2 + s 2 (33) 

and the kurtosis is 

^ = 3+5^(1 + ^)- (34) 

If fi2 = s 2 the nongaussian part of the kurtosis is reduced by a factor of four. 
This noise can either be instrumental or come from density perturbations at 
last scattering. Instrumental noise is uncorrelated, so that s 2 = 2a 2 if the 
noise has variance a 2 at each pixel. Gaussian noise from last scattering is 
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correlated, and some cold dark matter models predict the correlation function 
(Gott et al. 1990) 

C(0) = C(0)— V (35) 

1 + 2a 2 

for the temperature anisotropies (RV T) for scales 9 < 2°, where a = 8' is the 
coherence angle, and C(0) = a 2 is the variance of T. We can approximate the 
temperature difference between two neighbouring beams by Y = bT' and use 
the fact that if a gaussian RV T has the correlation function C(6), then its 
derivative T' has a gaussian distribution with variance — C"(0) (Vanmarcke 
1983). Since M bX (t) = M x {bt) , s 2 = -b 2 C"(0) = ^a 2 . A more careful 
analysis shows that this approximation is good for b < a. Fig. 8 shows the 
dependence of the kurtosis on noise from density perturbations at last scat- 
tering, and Fig. 9 for instrumental noise. For resolutions of a few arcminutes 
or below, noise from density perturbations at last scattering does not pose 
much of a problem because it is correlated, but instrumental noise can wash 
out the stringy signal if it is too large. 

As an example, we can use the above results to estimate the effect of 
cosmic string loops on the kurtosis of the CMB spatial gradient map. As 
discussed in section 2, loops contribute as gaussian noise via the Sachs- Wolfe 
integral. 

The Sachs- Wolfe integral due to cosmic string loops has been estimated 
in Brandenberger and Turok (1986) (see also Traschen, Turok and Bran- 
denberger (1986)). In that work it was shown that the largest loops at t\ s 
dominate the Sachs- Wolfe integral. The result for AT (6), the temperature 
difference between two points separated by an angle 9, scales as a 1//4 , where 
at is the radius of the largest loops at time t. Taking a = 1 and using 
G/z = 2 • 10~ 6 gave the value 

((^) 2 ) V2 ~5-10-W/^ £<0.5° (36) 

For 9 = 18", corresponding to the beam width of the VLA experiment, and 
with a = 10~ 2 we obtain a value of 1.5 ■ 10~ 7 . So the gaussian RV Y for the 
temperature gradients from loops has variance 

2 / ( AT X 2i 



s ={{-jr) )( d = b ) (37) 
From Fig. 9 it follows that the effect of loops on the kurtosis is negligible. 
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7 Discussion 



We have looked at a statistic which can distinguish strings from inflation and 
which can also be calculated analytically in a model for the string network, 
showing the dependence on string parameters. This is useful as a complemen- 
tary approach to numerical simulations which have to deal with singularities 
in the evolution equations of the strings leading to some uncertainty about 
the results. The statistic investigated here might not be the most sensitive 
one, we intend to look at others in the future. The genus curve (Gott et al. 
1990) for example, a topological statistic, makes use of the two-dimensional 
pattern of the anisotropies. A multifractal analysis (Pompilio 1993) is sensi- 
tive to higher moments of the temperature gradients. 
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9 Figure Captions 



Figure 1 : Simulated map of CMB temperature anisotropies from strings for 
a resolution of 80". The map has size 2.2° x 2.2° (100 x 100 pixels). There are 
M=10 strings per Hubble volume, and units are such that (3 = 47rGfi , ~f s v s = 1. 
The colour scale ranges from 5T/T = -15 to +15. 

Figure 2 : Model for the effect of a single string on the CMB temperature 
of a patch of sky of size 0x0 mapped out by circular beams. Strings located 
within the big square of size (0 + Oh) 2 can affect the temperature of the small 
patch. 

Figure 3 : Calculation of temperature difference of two adjacent beams 
due to one string. 

Figure 4 : Frequencies of temperature differences X (in units where 
P — 1) for b — 18" and a binsize of 0.2 for strings. The gaussian fit is chosen 
to have the same variance as the result for strings. 

Figure 5 : Kurtosis of the distribution of temperature differences be- 
tween adjacent beams expected for strings (with M = 10) as a function of 
beam size b. 

Figure 6 : Probability density functions of kurtosis of temperature gra- 
dients measured in the experiment by Fomalont et al. for an underlying 
stringy (for M = 10) and gaussian distribution. 

Figure 7 : Probability density of kurtosis measured in the experiment 
by Lasenby at 10' for strings (for M — 10 and M — 3) and an underlying 
gaussian distribution. 

Figure 8 : Dependence of kurtosis of temperature gradients on noise from 
density perturbations at last scattering giving rise to temperature anisotropies 
with variance a 2 at each pixel. Here and in the next figure the values M = 10, 
Gji = 2 x 10~ 6 and r y s v s = 0.15 are used. 
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Figure 9 : Dependence of kurtosis of temperature gradients on instru- 
mental noise (with variance a 2 at each pixel). 
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